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In the spacetime of n-dimensional static charged black hole we examine the mechanism by which 
the self-interacting scalar hair decay. It is turned out that the intermediate asymptotic behaviour 
of the self-interacting scalar field is determined by an oscillatory inverse power law. We confirm our 
results by numerical calculations. 
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■ (both analytically and numerically). It plays an important role in black hole's physics. Black hole radiates away 



The problem of the late-time behaviour of various fields in the spacetime of a collapsing body is well established 



everything that it can. This phenomenon happens regardless of details of the collapse or the structure and properties 
O i' of the collapsing body. The resultant black hole can be described only by few parameters such as mass, charge and 
angular momentum and due to the Wheeler's metamorphic dictum black holes have no hair (see for the vast amount 
^ ■ of references concerning this problem [1]). Therefore it is interesting to investigate how these hair loss proceed 
' dynamically. 

The neutral external perturbations were first studied in Ref. [2]. It was found that the late-time behavior is 
dominated by the factor ^-( 2i + 3 ) j f r each multipole moment I. On the other hand, the decay-rate along null infinity 
and along the future event horizon was governed by the power laws w~( i+2 ) and v~( l+3 \ where u and v were the 
outgoing Eddington-Finkelstein (ED) and ingoing ED coordinates [3]. The scalar perturbations on Reissner-Nordtrom 
(RN) background for the case when | Q |< M has the following dependence on time £~( 2 '+ 2 ) ; while for | Q \= M the 
late-time behavior at fixed r is governed by t~^ +2 ^ [4]. 

It turns out that a charged hair decayed slower than a neutral one [5]- [7], while the late-time tails in gravitational 
collapse of a self-interacting (SI) fields in the background of Schwarzschild solution was reported by Burko [8] and 
in RN solution at intermediate late-time was considered in Ref. [9]. At intermediate late-time for small mass m the 
decay was dominated by the oscillatory inverse power tails £ _ ('+ 3 / 2 ) sin(mi). This analytic prediction was verified 
at intermediate times, where mM < mt < l/(mM) 2 . In Ref. [10] the nearly extreme RN spacetime was considered 
and it was found analytically that the inverse power law behavior of the dominant asymptotic tail is of the form 
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t 5 / 6 sin(mi), independent of I. The asymptotic tail behaviour of SI scalar field was also studied in Schwarzschild 
spacetime [11]. The oscillatory tail of scalar field has the decay rate of i~ 5 / 6 at asymptotically late time. The power- 



law tails in the evolution of a charged massless scalar field around a fixed background of dilaton black hole was studied 
in Ref. [12], while the case of a self-interacting scalar field was elaborated in [13]. 

Nowadays it seems that that it is impossible to construct a consistent theory unifying gravity with other forces 
in Nature in four dimensions. The no-hair theorem for n-dimensional static black holes is quite well established 
[14]. So it will be not amiss to ask about the mechanism of decaying black hole hair in higher dimensional static 
black hole case. The evolution of massless scalar field in the n-dimensional Schwarzshild spacetime was determined 
in Ref. [15]. It was found that for odd dimensional spacetime the field decay had a power falloff like t~( 2l+n ~ 2 \ 
where n is the dimension of the spacetime. This tail was independent of the presence of the black hole. For even 
dimensions the late-time behaviour is also in the power law form but in this case it is due to the presence of black 
hole t-( 2( +3"-8) Gravitational perturbations of maximally symmetric black hole spacetime in higher dimensions 
were studied by Kodama et al. [16], while gravitational quasi- normal radiation of higher dimensional black holes was 
elaborated in [17]. 

In our work we shall consider and discuss the SI scalar field behaviour in the spacetime of n-dimensional static 
charged black hole. In Sec. II we gave the analytic arguments concerning the intermediate behavior of SI scalar field 
in the background of the considered black hole. Then, in Sec. Ill we treated the problem numerically and check our 
analytical considerations. We conclude our investigations in Sec. IV. 

II. MASSIVE SCALAR FIELDS IN N-DIMENSIONAL SPHERICALLY SYMMETRIC SPACETIME 

In our paper we shall investigate the evolution of SI (massive) scalar field tp in a fixed spacetime of a static electrically 
charged n-dimensional black hole. The wave equation for the field is given in the form as follows: 



where m is assumed to be real. 

The metric of the external gravitational field will be given by the static, spherically symmetric solution of equations 
of motion derived from the action written as 



where we denoted the generalized (n — 2)-gauge form F , Ml ... Alre _ 2 by the expression F( n -2) = dAr n _^\. Further, we 
assume that we have to do with the only one nontrivial electric component of the (n— 2)-gauge form as Aoi... n —2 = 4'{ x )- 
The metric of the spherically symmetric charged static n-dimcnsional black hole implies 
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where e?£l 2 



■n-2 



is the line clement of the unit S' 



'n-2 



sphere. Next, we define the tortoise coordinates y as 
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dy 



Thus, the metric (3) can be rewritten in the form as 



dr 



2M , Q 2 

r n-J I r 2n-ti 



(4) 



ds z = 1 - 
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- dt 2 + dy 2 
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(5) 



In the spherical background each of the multipole of perturbation field evolves separetly. Because of the fact that the 
scalar field is of the form 



(6) 



X 2 



where Y™ is a scalar spherical harmonics on the unit (n— 2)-sphere and m denotes a set of (n— 3) integers (mi . . . m„_ 3 ) 
satisfying I > m„_3 > . . . > m 2 >| mi |. In the process of this one has the following equations of motion for each 
multipole moment 



1p,tt ~ 1p,yy + Vl}) = 0. 



The potential V implies 



V = f 2 (r) 



T 2 



dr 
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where f 2 (r) = 1 

In order to analyze the time evolution of SI scalar field in the background of the considered black hole we shall use 
the spectral decomposition method [18]. The time evolution of SI scalar field may be written in the following form: 



G(y, y'; t)Mv' ', 0) + G t (y, y'; t)^(y', 0) 



V(»,t) - J dy 1 
for t > 0, where the Green's function G(y,y';t) is given by 

G(y,y';t) = 6(t)d(y-y'). 



dt 2 dy 2 



(9) 



(10) 



Our main task will be to find the black hole Green function so in the first step we reduce equation (10) to an ordinary 
differential equation. To do it one can use the Fourier transform [19] G(y, y 1 ; u) = J^_° dtG(y, y'; t)e luJt . This Fourier's 
transform is well defined for Im uo > 0, while the corresponding inverse transform yields 

G(y, y'; t) = — / dcu G(y, y'; u)e~^, (11) 

for some positive number e. The Fourier's component of the Green's function G(y,y';u)) can be written in terms of 
two linearly independent solutions for homogeneous equation as 

d 2 



dy 2 



+ cu z -V Vi = 0, i = 1,2, 



(12) 



The boundary conditions for ipi are described by purely ingoing waves crossing the outer horizon H + of the n- 



dimensional static charged black hole xpi ~ e 



as y — > — oo while ip2 should be damped expotentially at i + , namely 



v at y — > oo. 
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In order to find ipi we consider the wave Eq.(12) of SI scalar field and introduce an auxiliary variable £ in such a 
way that £ = f(r)ip. So in terms of £ relation (12) can be written as follows: 



d 2 i /(r), 



dr 2 f(r) 



f\r) p(r) 



n-2\ 1 d 



2 J j* 2 dv 



2 / 2 (r) ) + ^ + "-3) +m2 



(13) 



Further on, we expand Eq.(13) in power series of -j^ and neglecting terms of order Ol (-£=5) )■ Then, we 



arrive at the following expression: 

a 

dr 2 



where v = l-2+%. 



+ 



9 , AMuo 2 - 2Mm 2 v(v + 1) 
- 2 -- 2 + ^3 



If one further assumes that the observer and the initial data are in the region where r n < (Mm)' 1 an( ^ one 

M 
(Mm)' 

d 2 ^ 



interested in the intermediate asymptotic behavior of SI scalar field r n 3 <t< (h ¥ Yi , then we get 



dr 2 



2 2 v{v + l) 
lj — m k 



(14) 



shall be 



(15) 



As in four-dimensional case the scalar field perturbations on n-dimcnsional charged static black hole background does 
not depend on the spacetime parameters such as M and Q. The perturbations in question depend on the scalar field 
parameter (mass of the field) . 

The same procedure as described in Ref. [9] leads us to the solution of the relation (15) (we refer the readers to 
this work). Thus, in our case the intermediate asymptotic behaviour of the SI field at fixed radius has the form 



^~t-('+9-*), 



(16) 



while the intermediate behaviour of SI fields at the outer horizon H + is dominated by an oscillatory power law tails 
of the form as follows: 



ij; ~ v ^ + 2 2) sin(mt). 

In the next section we check our predictions numerically for various dimensions of the background spacetime. 
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III. NUMERICAL RESULTS 



We numerically analyzed Eq. (7) using method described in [3]. We transformed Eq. (7) into (u,v) coordinates 

4V>« + Vyj = 0, (18) 

and solved it on uniformly spaced grid using explicit difference scheme. As was pointed out previously the late time 
evolution of a massive field is independent of the form of the initial data. In order to perform our calculation we start 
with a Gaussian pulse of the form as follows: 

(v — v ) 2 \ , , 

(19) 



ip(u = 0, v) = A cxp — 
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Because of the linearity of the relation (7) one has freedom in choosing the value of the amplitude A. For our purpose 
we fix it as A = 1. The rest of the initial field profile parameters we take as vq = 50 and a = 2. 

We shall set the mass of SI scalar field equal to to — 0.01 and the mass and charge of the black hole respectively equal 
to M = 0.5, Q = 0.45. First, we shall study the evolution of ip on the future timelike infinity i + . In our calculations 
we approximate this situation by the field at fixed radius y — 50. The numerical results for / = and different 
spacetime dimensions n — 4, 5, 6 are shown in Fig. 1. Initially the evolution is determined by the prompt contribution 
and quasinormal ringing. However, then with the passage of time a definite oscillatory power-law fall off appears 
to be manifest. We obtained the power-law exponents —1.51, —2.03, and —2.53 for n = 4,5,6, respectively. These 
values are to be compared with the analytically predicted ones equal respectively to —1.5, —2.0, and —2.5. Thus, 
the agreement between numerical calculations and analytically predicted values is excellent. For the four-dimensional 
spacetime we get the same result as obtained in Ref. [9]. The period of the oscillations is T — tt /to ~ 314.5 ± 0.5 
to within 0.3% for all curves. Both values of power-law exponents and period of oscillation arc in perfect agreement 
with the analytical prediction. 

Then, the evolution of the SI scalar field on the black hole future horizon H + (approximated by tp{u — 10 4 , t)), as a 
function of v for different space dimensionality n — 4, 5, 6 was studied. Calculation parameters are: I = 0, M = 0.5, 
Q = 0.45, and to = 0.01. The power-law exponents and the period of the oscillations are the same as in Fig. 1. 
Next, we take into consideration the dependence of SI scalar field on the multiple index. We studied the evolution of 
the field %[> on the future timelike infinity i+ as a function of t for different multipoles I = 0, 1 and 2 in five-dimensional 
spacetime. The obtained power-law exponents are as follows: —2.03, and —3.03, and —4.04 for 1 = 0,1, 2. According 
to relation (16) these exponents are equal to —2.0, —3.0, and —4.0, in perfect agreement with numerical calculations. 
The period of the oscillations is T = ir/m ~ 314.5 to within 1% for all curves, in agreement with the predicted value. 
The results are depicted in Fig. 3. The same calculations were conducted for six-dimensional spacetime (Fig. 4), 
where the power-law exponents are —2.53, and —3.53, and —4.53 for I = 0,1,2, respectively. Due to Eq.(16) for 
six-dimensional spacetime they have the values: —2.5, —3.5, and —4.5. Those values are also in agreement with 
numerical calculations. The period of the oscillations is T = ir/m ~ 314.5 to within 3% (for the worst case of I = 2). 
We also studied numerically the late-time behaviour of SI on black hole in five-dimensional spacetime, for the field 
mass to = 0.05 (see Fig. 5). We investigated the behaviour on future timclikc infinity and on black hole future horizon. 
The period of oscillation was T = 63.0 to within 0.8%. The field's amplitude decays in agreement with the no-hair 
theorem due to the power-law fall off, contrary to the four-dimensional case where one has the decay rate slower than 
any power-law. In six-dimensional case ( Fig. 6) the behaviour of SI scalar fields is the same, with the period of 
oscillation equal to T = 63.0 to within 0.8%. The slope of the curve is equal to —2.5. 

We also investigated the late-time behaviour on future timelikc infinity in five-dimensional spacetime of static charged 
black hole for different masses of the scalar field Fig. 7. The decay rates have the form of the power-law fall-off with 
the slope equal to —2.0. The same studies were conducted in six-dimensional spacetime with the similar results Fig. 8, 
i.e., for various masses of SI scalar fields we obtained power-law fall off with the slope of the curve equal to —2.5. 
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IV. CONCLUSIONS 



We have elaborated analytically the intermediate behaviour of SI scalar fields in the background of static charged 
n-dimensional black hole. It turned out that the intermediate asymptotic behaviour did not depend on the spacctimc 
parameters as M, Q but only on the mass of SI scalar field. In other words, considering the SI scalar field perturbations 
one can neglect the backscattcring from the asymptotically far regions at intermediate times. We check our analytical 
calculations by numerical ones and obtained excellent agreement with analytically predicted values. Numerical studies 
of the late-time behaviour of SI scalar fields on the future timelike infinity and on the future black hole horizon reveal 
the fact of the power-law decay (contrary to the four-dimensional case where one has to do with the slower than 
any power-law decay). The same type of behaviour was obtained for various masses of SI scalar fields in n = 5,6 
spacetime. This type of behaviour should be checked analytically, but due to the tremendous difficulties in solving 
differential equations is almost intractable. We hope to return to this problem elsewhere. 
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FIG. 1. Evolution of the field \tp\ on the future timelike infinity i + (approximated by ip(y = 50, £)) as a function of t for 
different space dimensionality n — 4, 5, 6 (curves from top to bottom, respectively). Calculation parameters are: I = 0, M = 0.5, 
Q = 0.45, and m = 0.01. The powerdaw exponents are —1.51, —2.03, and —2.53 for n — 4,5, 6, respectively. The period of the 
oscillations is T — ir/m ~ 314.5 ± 0.5 to within 0.3% for all curves. 
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FIG. 2. Evolution of the field \ip\ on the black hole future horizon H+ (approximated by tp(u = 10 4 , t)) as a function of v 
for different spacetime dimensionality n = 4, 5, 6 (curves from top to bottom, respectively). Calculation parameters are: I — 0, 
M — 0.5, Q — 0.45, and m = 0.01. The power-law exponents and the period of the oscillations are the same as in Fig. 1. 
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FIG. 3. Evolution of the field \tp\ on the future timelike infinity i + (approximated by ip(y = 50, t)) as a function of t for 
different multipoles I — 0, 1 and 2 (curves from top to bottom, respectively) and n = 5. The power-law exponents are —2.03, 
and —3.03, and —4.04 for I = 0,1, 2, respectively. The period of the oscillations is T = ir/m ~ 314.5 to within 1% for all curves. 
The rest of the parameters are the same as in Fig. 1. 
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FIG. 4. Evolution of the field \tp\ on the future timelike infinity i + (approximated by ip(y = 50, t)) as a function of t for 
different multipoles I — 0, 1 and 2 (curves from top to bottom, respectively) and n = 6. The power-law exponents are —2.53, 
and —3.53, and —4.53 for I = 0, 1, 2, respectively. The period of the oscillations is T = ir/m ~ 314.5 to within 3% (for the 
worst case of I = 2). The initial data are those of Fig. 1. 
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FIG. 5. Late-time behaviour of the field \ip\ on black hole for n = 5. Only maxima of the oscillations are shown. The mass 
of the field is m = 0.05. The upper curve represents the field on future timelike infinity i + (if>{y = 50)) as a function of time t. 
Bottom curve is the field on the black hole future horizon H + (ip(u = 2 x 10 4 ) as a function of v. The period of the oscillations 
is T = 7r/m ~ 63.0 to within 0.8%. The thin dashed lines have slopes equal to —2.0. The rest of the parameters are the same 
as in Fig. 1. 
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FIG. 6. Late time behaviour of the field on black hole for n = 6. Only maxima of the oscillations are shown. The mass 
of the field is m = 0.05. The upper curve represents the field on future timelike infinity i + (V>(j/ = 50)) as a function of time t. 
Bottom curve is the field on the black hole future horizon H + (tp(u = 2 x 10 4 ) as a function of v. The period of the oscillations 
is T = -iv/m ~ 63.0 to within 0.8%. The thin dashed lines have slopes equal to —2.5. The rest of the parameters are the same 
as in Fig. 1. 
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FIG. 7. Late-time behaviour of the field \ip\ max at future timelike infinity on black hole for n = 5 and for different masses 
m. The dashed line has slope equal to —2.0. The initial data are those of Fig. 1. 
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